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We investigate the nature of the near-threshold enhancement in the pp invariant mass spectrum 
of the reaction J/ty— >~/pp reported recently by the BES Collaboration. Using the Julich NN model 
we show that the mass dependence of the pp spectrum close to the threshold can be reproduced by 
the S-wave pp final state interaction in the isospin 1=1 state within the Watson-Migdal approach. 
However, because of our poor knowledge of the NN interaction near threshold and of the J/ty-^ypp 
reaction mechanism and in view of the controversal situation in the decay J/9— >ir°pp, where no 
obvious signs of a pp final state interaction are seen, explanations other than final state interactions 
cannot be ruled out at the present stage. 
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I. INTRODUCTION 

Recently the BES Collaboration [1] reported a near- 
threshold enhancement in the proton-antiproton (pp) in- 
variant mass spectrum, observed in the J/ty—*jpp decay. 
Signs for a low mass pp enhancement had been already 
seen earlier by the Belle Collaboration in their study 
of the B + ^K + pp decay [2] as well as in the reaction 
B°^D°pp [3]. But because of the large statistical un- 
certainties of the Belle data it was difficult to draw con- 
crete, quantitative conclusions about the extent of the 
near threshold pp enhancement. The new data by Bai 
et al. [1], however, are of rather high statistical accuracy 
and therefore provide very precise information about the 
magnitude and the energy dependence of the pp mass 
spectrum very close to threshold. 

The BES Collaboration [1] htted their pp invari- 
ant mass spectrum below 1.95 GeV by a Breit-Wigner 
resonance function. Assuming that the pp system 
is in an S'-wave resulted in a resonance mass of 
M=1859t?^ 25 MeV and a total width of T<30 MeV. 
A comparable fit to the data could be achieved with a 
P-wave Breit-Wigner function with A/=1876±0.9 MeV 
and r=4.6±1.8 MeV. 

The proximity of these resonance masses to the pp re- 
action threshold (which is at 1876.54 MeV) nourished 
speculations that the observed strong enhancement could 
be a signal of an NN bound state. While theoretical 
considerations of such NN bound states (or of baryonia, 
in general) abound in the literature [4-9] there is so far 
hardly any undisputed experimental information on the 
existence of such states [10-12]. Thus, the supposition 
that one has found here independent and possibly even 
more convincing evidence in support of NN bound states 
is certainly appealing. 

An alternative explanation put forward by the BES 
collaboration invokes similarities of the observed en- 
hancement in the pp mass spectrum near threshold with 



the strong energy dependence of the electromagnetic 
form factor of the proton around y/s~2m p , in the time- 
like region, as determined in the reaction pp—>e + e~ [13]. 
In the latter case it was argued that the sharp struc- 
ture seen in the experiment could be caused by a nar- 
row, near-threshold vector-meson (J PC =1 ) resonance 
[14] with M=1870±10 MeV and T=10±5 MeV. (See also 
Refs. [15, 16] for a pertinent discussion.) One should 
keep in mind, however, that the quantum numbers of 
the J/\E' particle (J PC =1~~) would restrict such a reso- 
nance to occur in a pseudoscalar (0~~) or scalar (0 H ) 
state - should it be indeed responsible for the enhance- 
ment of the near-threshold pp mass spectrum in the decay 
J / ty^jpp. 

An entirely different and much more conventional in- 
terpretation of the observed enhancement was suggested 
in several recent works [17-19]. These authors argue 
that the enhancement is primarily due to the final state 
interaction (FSI) between the produced proton and an- 
tiproton. Specifically, it was shown within the scattering 
length approximation [17] that a calculation with a com- 
plex S-w&ve scattering length extracted from an effective- 
range analysis of pp scattering data can reproduce the 
shape of the pp mass distribution close to the threshold. 

In the present paper we analyse the near threshold 
enhancement in the pp invariant mass spectrum reported 
by the BES Collaboration utilizing a realistic model of 
the JVJV interaction [21]. The elastic part of this model is 
the G-parity transform of the Bonn meson-exchange JVJV 
potential, supplemented by a phenomenological complex 
potential to account for NN annihilation. 

As just mentioned, the investigations of Kerbikov et al. 
[17] as well as those of Bugg [18] rely on the scattering 
length approximation. But it remains unclear over which 
energy range this approximation can provide a reliable 
representation of the energy dependence of the pp ampli- 
tude. The enhancement seen in the BES data extends 
up to invariant masses of M(pp)«2 GeV, which corre- 
sponds to center-of-mass energies of around 120 MeV in 
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the pp system. It is obvious that the simple scattering 
length approximation cannot be valid over such a large 
energy region. Using the pp amplitude of our interac- 
tion model (which describes the available NN scattering 
data up to ccntcr-of-mass energies of 150 MeV [21]) we 
can examine to what extent the scattering length approx- 
imation can indeed reproduce the energy dependence of 
the pp scattering amplitude. Moreover, we can compare 
the energy dependence induced by the full pp amplitude 
with the BES data over the whole energy range where 
the enhancement was observed. 

A microscopic NN model has a further advantage. It 
yields predictions for all possible spin- and isospin chan- 
nels. Thus, we can investigate whether the energy de- 
pendence of the 1 S'o and 3 Po NN scattering amplitudes 
in the 1=0 and 7=1 isospin channels is compatible with 
the BES data. (We use the nomenclature ( 2 '+i)(2S+i) Lj 
but omit the isospin index when referring to both isospin 
channels or to isospin averaged results.) The analysis of 
Kerbikov et al. utilizes scattering lengths obtained from 
a spin-averaged effective-range fit to the low-energy NN 
data. Those values might be strongly influenced or even 
dominated by the 3 Si NN partial wave, a channel which 
cannot contribute to the FSI in the J/*— >7pp decay be- 
cause of charge-conjugation invariance. 

Let us finally mention that the work of Zou and Chiang 
[19] does not employ the scattering length approximation 
but uses the if-matrix formalism. However, their inter- 
action consists only of one-pion exchange and is there- 
fore not realistic because it does not take into account 
the most striking feature of low-energy NN scattering, 
namely annihilation. Recall that near the pp threshold 
the annihilation cross section is twice as large as the elas- 
tic pp cross section [21]. 

The paper is structured in the following way: In 
Sect. II we briefly review the NN potential model that is 
used in the present analysis and examine the reliability of 
the scattering length approximation on the NN ampli- 
tude for the partial waves relevant to the analysis of the 
BES data. The possible influence of the Coulomb inter- 
action is discussed as well. In Sect. Ill we provide details 
of our calculation of the near-threshold pp mass spectrum 
for the reaction and we compare our results with the mea- 
surement of the BES collaboration. Sect. IV is devoted 
to a discussion of possible signals of NN bound states 
or (sub pp threshold) meson resonances in the pp mass 
spectrum. The paper ends with concluding remarks. 



II. PROTON- ANTIPROTON SCATTERING AT 
LOW ENERGIES 

In the present investigation we use one of the NN mod- 
els developed by the Jiilich group [20-22]. Specifically, 
we use the Bonn OBE (one-boson-exchange) model in- 
troduced in Ref. [20] whose results are also discussed in 
Ref. [21]. In the latter work the model is called A(OBE) 
and we also adopt this name in the present paper. The 
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FIG. 1: The pp elastic cross section as a function of the invari- 
ant collision energy. The lines show the result of the Jiilich 
model A(OBE) [21]. The dashed line is the contribution from 
the 1 So partial wave alone, dotted - from all s waves, dashed- 
dotted from p waves and solid lines is full calculations. The 
circles are experimental data taken from Ref. [47] . 



model A(BOX) is based on the OBE version of the full 
Bonn potential derived in time-ordered perturbation the- 
ory, i.e. OBEPT in Ref. [23]. The G-parity transform of 
this interaction model constitutes the elastic part of the 
NN potential. A phenomcnological spin-, isospin- and 
energy-independent complex potential of Gaussian form 
is added to account for the NN annihilation. With only 
three free parameters (the range and the strength of the 
real and imaginary parts of the annihilation potential) 
a good overall description of the low- and intermediate 
energy NN data was achieved. Results for the total and 
the integrated elastic and charge- exchange cross sections 
as well as angular dependent observables can be found in 
Refs. [20, 21]. Here in Fig. 1 we show only the elastic pp 
cross section. 

We also utilized the most complete NN model of the 
Jiilich Group, model D published in Ref. [22]. The elastic 
part of this interaction model is derived from the G-parity 
transform of the full Bonn NN potential and the anni- 
hilation is described in part microscopically by NN^ 2 
meson decay channels - see Ref. [22] for details. But since 
the results turned out to be qualitatively rather similar to 
the ones obtained with A(OBE) we refrain from showing 
them here. 

The differential cross sections for pp scattering [24-26] 
already indicate a substantial contribution from higher 
partial waves at around 10 MeV above the pp threshold. 
This feature is also reflected in the predictions of the NN 
model. The dashed line in Fig. 1 shows the contribution 
from the x Sq partial wave. The dotted and dash-dotted 
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FIG. 2: The pp scattering amplitudes for the 1 So and 3 Si par- 
tial waves as a function of the center-of-mass energy. The solid 
lines show the results of the Jiilich model while the dashed 
lines indicate the scattering length approximation given by 
Eq.(l). 



lines illustrate the total contributions from s and p waves, 
respectively. 

Since the scattering length approximation to the NN 
amplitude was used in two of the analyses of the BES 
data [17, 18] we want to investigate the validity of this 
approximation. In the following wc ignore the proton- 
neutron mass difference and also the Coulomb interaction 
in the pp system, in order to simplify the discussion. But 
we will come back to these issues later. 

In the scattering length approximation the S-wetve pp 
scattering amplitude T is given by 



T = 



iaq p 



(1) 



where the scattering length a is a complex number be- 
cause of inelastic channels (annihilation into multime- 
son states) that are open already at the pp thresh- 
old. The proton momentum in the cm. system is 



4m2/2, where m p is the proton mass. 

Results for \T\ 2 for the 1 S'o partial wave are presented 
in the upper part of Fig. 2. The solid lines are the re- 
sult for the full amplitude while the dashed lines are 



based on the scattering length approximation given by 
Eq.(l). Note that the scattering lengths predicted by 
the NN model, which we use for the partial wave 
are do=(— 0.18— il.18) fm and di=(1.13— i0.61) fm for 
the isospin 1=0 and 1=1 channels, respectively. It is ev- 
ident that the scattering length approximation does not 
reproduce the energy dependence of the scattering am- 
plitude that well. For the 7=1 channel the difference at 
an excess energy of 50 MeV amounts as much as 50 %. 
The difference is even more pronounced for the 1=0 chan- 
nel, where we already observe large deviations from the 
full result at rather low energies. This strong failure of 
the scattering length approximation is due to the much 
smaller scattering length predicted by our model for the 
1=0 partial wave. 

Results for the 3 Si partial wave are shown in the 
lower part of Fig. 2. Here the scattering lengths pre- 
dicted by the NN model are ao=(1.16— i0.82) fm and 
ai = (0.75— i0.84) fm for the 1=0 and 1=1 channels, re- 
spectively. As pointed out above, this partial wave can- 
not contribute to the reaction J/^^ypp. But it is still 
interesting to see how well the scattering length approxi- 
mation works in this case. Obviously we also see similar 
shortcomings here. Thus, it is clear from Fig. 2 that the 
use of the scattering length approximation allows only a 
very rough qualitative estimate of the effects of FSI but it 
is definitely not reliable for a more quantitative analysis 
of the BES data. 

Nevertheless, one has to realize that the main uncer- 
tainty in estimating pp FSI effects does not come from 
the scattering length approximation but from our poor 
knowledge of the pp 1 Sq amplitudes near threshold and 
of the J/^— >7pp reaction mechanism. For example, the 
scattering lengths employed by Kerbikov et al. are spin- 
averaged values. Since the 3 Si partial wave contributes 
with a weighting factor 3 to the pp cross sections (and 
there are no spin-dependent observables at low energies 
that would allow one to disentangle the spin-dependence) 
it is obvious that their value should correspond predom- 
inantly to the 3 5*i amplitude. Thus, it is questionable 
whether it should be used for analyzing the BES data at 
all because the contribution of the 3 <Si partial wave to the 
decay J/^— ^pp is forbidden by charge-conjugation in- 
variance. But even the availability of genuine 1 Sq ampli- 
tudes (which are necessarily model dependent) does not 
solve the problem. The reaction J/^=>jpp can have any 
isospin combination in the final pp state. A dominant 
reaction mechanism involving an intermediate isoscalar 
meson resonance (i], r] c , ■■•), cf. Fig. 3c, would yield a 
pure 1=0 final pp state, whereas an intermediate isovec- 
tor meson resonance (jr, 7r(1800), ...) leads to pure 1=1. 
On the other hand, mechanisms like the ones depicted in 
Fig. 3a, b generate the usual equal weighting of the two 
isospin amplitudes. 

In this context let us mention that the spin-averaged 
pp scattering length predicted by the model A (BOX), 
a5=0.84— iO. 85 fm, is in rough agreement with the 
value extracted from the level shifts of antiproton-proton 
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a large scattering length similar to the one in the 1=1 
channel, whereas modifications of the short-range con- 
tributions to the elastic part of the TV TV potential had 
hardly any qualitative influence on the scattering lengths. 
Thus, it would be interesting to obtain experimental con- 
straints on the 1 So scattering lengths. Corresponding 
measurements could be performed at the future GSI fa- 
cility FLAIR [30] , where it is possible to have a polarized 
antiproton beam [31], as required for inferring the spin- 
triplet and spin-singlet amplitudes from the data. 



7,7T 




M 



J/M> 
c) 



J/* 
d) 



FIG. 3: Some possible reaction mechanisms for the de- 
cays J/ty— >7pp and J/vp— *Tv°pp. M indicates intermediate 
mesonic states and N (N) intermediate nucleon (antinucleon) 
resonances. 



atoms, cf. Table 8 in Rcf. [27]. On the other hand, this 
averaged value differs significantly from the scattering 
length in the 11 So partial wave - which is one of the possi- 
ble final states in the decay J/ , I>— >7pp. At first glance one 
would tend to believe that the unusually small real part 
of ao is basically accidental and reflects the fact that the 
1 S'o state enters with smaller statistical weight into the 
pp cross section and, therefore, is much less constrained 
by the pp data. However, exploratory calculations with 
the other TV TV models of the Jiilich group revealed that 
the real part of the 11 So scattering length is always small 
(and often of a different sign from the other S waves). 
For example, the more involved NN model D published 
in Ref. [22], yields a = - 0.25-il.Ol fm. Moreover, the 
NN models of Dover- Richard also and Kohno-Weise pre- 
dict rather small values for the real part of the 11 Sb scat- 
tering length, as is evident from the results presented in 
Ref. [28]. The Paris NN potential, on the other hand, 
predicts the real parts of the 1 S'o scattering lengths to 
be of similar magnitude - cf . Table IV in Ref. [29] . But 
Fig. 11 of that paper suggests that the central potential 
in the J = 0, S = channel has been drastically mod- 
ified inside 1 fm as compared to what follows from the 
G-parity transformation of the Paris NN potential. In 
the Jiilich TV TV models no such changes are introduced 
when deriving the elastic part of the TV TV interaction via 
G-parity! 

A closer examination of our NN models disclosed 
that here the smallness of Re an is definitely caused by 
the one-pion exchange contribution, i.e. by long-range 
physics. Switching off this contribution always led to 



III. J/* DECAY RATE AND FSI EFFECTS 

The J/^— >7£>p decay rate is given as [32] 
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x A 1/2 (M 2 , m 2 , m 2 ) dMdfl p dQ 7 , 
where the function A is defined by 

(x — y — z) 2 — 4yz 



X(x,y,z) 
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(2) 



(3) 



M is the invariant mass of the pp system, fi p is the proton 
angle in that system, while f2 7 is the photon angle in the 
J/ty rest frame. After averaging over the spin states and 
integrating over the angles, the differential decay rate is 



dT 
dM 



(m 2 



M 2 



Tkf 2 



4m 2 



(4) 



where A is the total J/^^jpp reaction amplitude. Note 
that A is dimensionless. The differential rate dT/dM, in- 
tegrated over the range from 2m p to mjm, yields the par- 
tial J/^^pp decay width, r=(3.3±0.9)lCT 2 keV [47]. 

The solid line in Fig. 4 is obtained with Eq. (4) by 
using a constant |^4| 2 . This result corresponds to the so- 
called phase space distribution. The actual value of |^4| 2 
was adjusted to the data at M>2.8 GeV, i.e. around 
the r] c resonance region. The circles in Fig. 4 show data 
for the J/^—^-fpp decay published by the BES Collabo- 
ration [1]. Evidently, close to the pp threshold the data 
deviate substantially from the phase space distribution. 

Since close to the threshold the phase space factor in- 
troduces a strong but trivial energy dependence of the 
mass distribution it is convenient to divide the experi- 
mental pp spectrum by the kinematical factors that ap- 
pear on the right side of Eq.(4) and extract the invariant 
amplitude \A\ 2 from the data. Corresponding results for 
the BES data [1] are shown by the squares in Fig. 5. It 
is obvious that the experimental reaction amplitude ex- 
hibits a very strong mass dependence near the pp thresh- 
old. In recent investigations [17-19] this feature has been 
attributed to a strong final state interaction (FSI) be- 
tween the outgoing protons and antiprotons. Indeed, an 
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FIG. 4: The pp mass spectrum from the decay J/*S>—>-ypp. 
The circles show experimental results of the BES Collabora- 
tion [1], while the solid line is the spectrum obtained from 
Eq.(4) by assuming a constant reaction amplitude A. 



enhancement and/or a strong energy dependence of near- 
threshold cross sections has been seen in many other re- 
actions involving hadrons and is commonly seen as to be 
caused by FSI effects, cf. [33-39]. Thus, it is plausi- 
ble and it even has to be expected that FSI effects play 
a role in the decay J/'F— >"/pp as well. The interesting 
issue is primarily whether the proton- antiproton inter- 
action is really strong enough to induce such a sizeable 
enhancement in the invariant mass spectrum. 

A very simple and therefore also very popular treat- 
ment of FSI effects is due to Watson [40] and Migdal [41]. 
These authors suggested that the reaction amplitude for 
a production and/or decay reaction which is of short- 
ranged nature can be factorized in terms of an elementary 
production amplitude A and the pp scattering amplitude 
T of the particles in the final state, 



Aprod ~ NAq ■ T, 



where Ao is given by diagrams like those shown in Fig. 
3 and N is a normalization factor. 

If the production mechanism is of short range then 
the amplitude A depends only very weakly on the en- 
ergy and the near-threshold energy dependence of the 
reaction amplitude is driven primarily by the scattering 
amplitude, T, of the outgoing particles. This means that 
in the case of the reaction J/^^jpp the near threshold 
mass dependence of the pp spectrum should be domi- 
nated by the energy dependence of the pp scattering am- 
plitude. In the analyses of the BES data by Kerbikov 
et al. [17] and by Bugg [18] the above treatment of FSI 
effects was adopted. 



However, the prescription of Watson-Migdal is only 
valid for interactions that yield a rather large scattering 
length, like the 1 Sq NN partial wave where the scattering 
length a is in the order of 20 fm, and even then only for 
a relatively small energy range, as was pointed out in 
several recent papers [42-44]. Therefore, in the present 
case, where the scattering lengths are in the order of 
1 fm one should be cautious with the interpretation of 
results obtained from applying Eq. (5) . Rather one should 
start from the more general expression for the reaction 
amplitude, 



A 



prod 



A a + A a G pp T, 



(6) 



which corresponds to a distorted wave born approxima- 
tion. If one assumes that the production amplitude, 
Aprod, has only a very weak energy- and momentum de- 
pendence it can be factorized and one obtains 



A prod ~A [1 + G™T] = A,*(-)*(0), 



(7) 



where *g p ^*(r) is a suitably normalized NN continuum 

wave function and <Fq p '*(0) is nothing else but the inverse 

of the Jost function J, i.e. ^^(O)^"^-^) [45]. 
Eq.(6) itself can be cast into the form [43] 

A.„ r „d = A [l + (c - iq p )T], (8) 



where c is, in general, a complex number that represents 
the principal-value integral over the half-off-shell exten- 
sion of the production amplitude Aq and the pp scatter- 
ing amplitude that appears on the very right hand side 
of Eq.(6). Obviously, Eq.(8) is formally equivalent to the 
Watson-Migdal prescription of Eq.(5) if \c\ is large. In 
case of the NN 1 Sq partial wave |c| is indeed large, as 
has been demonstrated in Ref. [44]. 

In considering the S*-wave interaction in the final pp 
system one should account for the centrifugal barrier be- 
tween the photon and the pp state given by the factor 



M 2 



2m 



(9) 



(5) with lj=l being the orbital momentum between the pho- 
ton and the pp system. Around the threshold, M=2m p , 
the centrifugal correction depends only weakly on the 
invariant mass of the pp system and basically does not 
modify the mass dependence of the FSI. Within the range 
2m p <M<2.1 GeV the factor C\ varies between 0.23 and 
0.27, which is not large enough to counterbalance the 
contribution from the strong S'-wave FSI over the same 
invariant-mass range. 

Note also that the Watson-Migdal formula (5) has to 
be modified when used for pp FSI effects in a P-wave. 
Since in the production reaction only the final momentum 
is on-shell one has to divide the on-shell T matrix by 
the factor q p in order to impose the correct threshold 
behavior of the production amplitude. 
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IV. RESULTS 

Let us first discuss calculations based on the Watson- 
Migdal approach given by Eq. (5). The solid lines in Fig. 
5 show the pp invariant scattering amplitudes squared 
for the 1 So and 3 Pq partial waves and the 1=0 and 1=1 
channels. We consider the isospin channels separately 
because, as mentioned above, the actual isospin mixture 
in the final pp system depends on the reaction mecha- 
nism and is not known. Note that all squared pp scat- 
tering amplitudes \T\ 2 were normalized to the BES data 
at the invariant mass M(pp)—2m p =50 MeV by multiply- 
ing them with a suitable constant. The results indicate 
that the mass dependence of the BES data can indeed be 
described with FSI effects induced by the scattering 
amplitude in the 1=1 isospin channel. The 1=0 chan- 
nel leads to a stronger energy dependence which is not in 
agreement with the BES data. We can also exclude dom- 
inant FSI effects from the 3 Po partial waves. Here the 
different threshold behaviour due to the P-wave nature 
cannot be brought in line with the data points very close 
to threshold. It should be clear, of course, that a suitable 
combination of several partial waves might as well repro- 
duce the experimental results on the pp invariant mass 
spectrum. 

Note that we do not include the Coulomb interaction in 
our model calculation and we also ignore the difference 
in the pp and nn thresholds. Judging from the results 
shown by Kerbikov et al. [17] their influence is noticeable 
only for excess energies below say 5 MeV. Accordingly 
we do not consider the lowest data point of the BES 
experiment, which lies in the first (5 MeV) energy bin, 
in our discussion. 

Our results support the conjecture of Kerbikov [17] 
and Bugg [18] that the enhancement seen in the near- 
threshold pp invariant mass spectrum of the decay 
J/^^jpp could be primarily due to FSI effects in the 
pp channel. But we should also consider the pp invari- 
ant mass spectrum of the decay J/^=nT pp which was 
presented by the BES collaboration in the same paper 
[1]. Here the near-threshold mass distribution does not 
show any enhancement as compared to the phase-space 
distribution - though one would likewise expect strong 
FSI effects in the pp channel. In this reaction isospin is 
preserved and the possible partial waves in the final pp 
state, 33 Si and 31 Pi, differ from those available in the 
decay J/^~ >^pp. But the 33 Si partial wave amplitude 
of the Jiilich NN model, shown in the lower part of Fig. 
2, leads to FSI effects that are comparable to those of 
the 31 So state - and in contradiction with the experi- 
mental ir°pp mass distribution - as illustrated in Fig. 6. 
We should mention in this context that the issue of the 
TT°pp mass distribution is not discussed in the work of 
Bugg, while Kerbikov argues that in the effective range 
parameterization that he employs \Re ai \ is much smaller 
than \Re ao\ and therefore there should be much smaller 
FSI effects in the 1=1 channel. However, our experience 
is that pp amplitudes with a small \Re a\ often yield an 
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FIG. 5: Invariant J/ 1 J>— >7pp amplitude \A\ 2 as a function of 
the pp mass. The squares represent the experimental val- 
ues of | A | 2 extracted from the BES data[l] via Eq.(4). The 
curves are the scattering amplitude squared (|T| 2 ) predicted 
by the AW model A(OBE) for the and 3 P partial waves 
and the 1=0 (solid) and 1=1 (dashed) channels, respectively. 
Note that the latter results have been normalized to |^4| 2 at 
M(pp)-2mp=50 MeV. 



even stronger energy dependence and therefore larger FSI 
effects - see the case of the 11 So partial wave in Fig. 2. 

In our opinion these seemingly contradictory results 
can only be reconciled if we assume that the treatment of 
FSI effects by means of Eq.(5) is oversimplified and one 
should use Eq.(8) instead. Results based on the latter 
equation are shown by the solid lines in Fig. 6 where 
we set c= — 0.1. It might be possible that a fine tuning 
of c could allow one to better reproduce the BES data 
for the TT°pp channels, though in view of the rather large 
statistical variations in the data we refrain from doing so. 
In any case it is obvious that the more refined treatment 
of the FSI effects based on Eq. (8) leads, in general, to 
a weaker energy dependence of the pp mass spectrum. 



V. DISCUSSION 

It is important to note that the near-threshold con- 
tribution to the J /^=f^pp decay rate is relatively large. 
The contribution to J/^^jpp up to invariant pp masses 
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FIG. 6: (a) The pp mass spectrum from the decay J/^— >7r°pp. 
The circles show experimental results of the BES Collabora- 
tion [1]. (b) Invariant J /^P— >7r°pP amplitude \A\ 2 as a func- 
tion of the pp mass. The circles represent the experimental 
values of \A\ 2 extracted from the BES data via Eq. (4). The 
curves are corresponding calculations using the T matrix pre- 
dicted by the NN model A(OBE) for the 33 Si partial wave. 
The dashed line is the results from the Watson-Migdal pre- 
scription Eq. (5), while the solid line is based on Eq. (8) with 
c— — 0.1. Both curves are normalized arbitrarily. 



of M(pp)m 2 GcV is roughly 5 times larger than the rate 
for the reaction J/ty— >7rj c followed by the rj c ^pp decay. 
Taking into account that the latter rate was recently pub- 
lished [46] as BR=1.9T0 -5 we can estimate that 

BRiJ/^^So) x BRCSo^pp) 9.5 • 10" 5 , (10) 

where 1 Sq indicates the near-threshold contribution of 
this partial wave - which is indeed a large fraction 
of the total J/vI/— >7pp branching rate that amounts to 
(3.8±1.0)-10" 4 . 

This suggests that in this energy region the reaction 
J/^^jpp could be indeed dominated by a resonance. 
Such a resonance should lie below the pp threshold but, 
unlike the resonance state which emerged from the Breit- 
Wigner fit of the BES collaboration [1], its mass could be 
significantly below the pp threshold and it could have a 
large width, like the 7r(1800) or 77(1760) resonances that 
are listed in the Particle Physics Booklet [47] . We demon- 
strate the scenario in Fig. 7. For illustration we consider 
the 7r(1800) meson with mass and width 1801 and 210 
MeV, respectively. The corresponding resonance ampli- 
tude squared is shown by the dashed line in Fig. 7. Its 



energy dependence is, of course, too weak and does not 
agree with the near-threshold pp mass spectrum found by 
the BES collaboration. However, when the 7r(1800) res- 
onance decays into a proton and antiproton, a final state 
interaction should occur. The solid line in Fig. 7 shows 
the total reaction amplitude squared, taken now as the 
product of the resonance amplitude and the pp scattering 
amplitude in the 31 So partial wave as is given by Eq.(5). 
The corresponding result is roughly in line with the BES 
data. 

Though this line of reasoning is certainly more a plau- 
sibility argument, rather than a solid calculation, we be- 
lieve that it makes clear that one has to be rather cau- 
tious when trying to extract resonance properties by fit- 
ting the near threshold pp mass spectrum with a reso- 
nance amplitude, because any final state interaction will 
necessarily and substantially distort the production am- 
plitude. 

A different and alternative scenario consists in assum- 
ing the formation of many multimcsonic intermediate 
states in the J/^f decay that couple strongly to the pp 
system where again a strong FSI occurs. The branch- 
ing rates of the J/<F radiative decays into the ir + ir~2ir°, 
pp, 2ir + 2ir~ and ujlo channels are larger than 10~ 3 and 
the transition of those mesonic states into the pp final 
state could be sufficiently strong [21, 22]. Thus, such a 
scenario cannot be excluded by the available data. 

Finally, let us come back to the other proposed ex- 
planation for the enhancement in the pp mass spectrum, 
namely near-threshold NN bound states. Clearly a de- 
scription of the experimental mass spectrum in terms of 




2 2.1 
M(pp) (GeV) 

FIG. 7: Invariant J/*— >7pp amplitude \A\ 2 as a function of 
the pp mass. The squares represent the experimental values of 
\A\ 2 extracted from the BES data[l] via Eq.(4). The dashed 
line shows the 7r(1800) resonance amplitude squared, while 
the solid line indicates the square of the product of resonance 
and pp 1 So scattering amplitude given by Eq.(5). The arrow 
shows the pp threshold. 



8 



pp FSI effects does not contradict the existence of such 
states. Indeed, in case of the NN interaction the strong 
FSI effects are interconnected with the existence of a 
near-threshold bound state (deuteron) or anti bound- 
state in the corresponding 3 Si and 1 Sq partial waves. 
However, the AW model that we used in the present 
study does not lead to any near-threshold bound states. 
In fact, we found only two bound states for the model 
A(OBE) within 100 MeV from the threshold, namely at 
E = -104-i413 MeV in the n S partial wave and at E 
= -24.2-H07 MeV in the 13 P partial wave. Obviously 
both bound states lie rather far away from the real axis 
and therefore should have practically no influence on the 
pp scattering amplitude near threshold. 

VI. CONCLUSION 

We have investigated suggested explanations for the 
near-threshold enhancement in the pp invariant mass 
spectrum of the J/*— >7pp decay, reported recently by 
the BES Collaboration. In particular, we showed that 
the near-threshold enhancement in the pp mass spec- 
trum can, in principle, be understood in terms of a final 
state interaction in the outgoing proton-antiproton sys- 
tem. Within our model calculation it can be described 
with the scattering amplitude in the So pp partial wave 
but disagrees with the energy dependence that follows 
from the n So and both (I = 0, 1) 3 Po amplitudes. 

We showed that the scattering length approximation 
cannot be used for a more quantitative evaluation of the 
energy dependence of BES data. In general, it does not 
reproduce the energy dependence of the scattering ampli- 
tude reliably enough within the pp mass range required. 
At the same time, one has to concede that the lack of 
knowledge of the reaction mechanism for J/^^jpp and 
the insufficient information on the pp interaction near 
threshold precludes any more quantitative conclusions, 
even when microscopic models are used - as in the case 
of the present study. 

We argued also that the use of the simple Watson- 
Migdal prescription for treating FSI effects has to be con- 
sidered with caution for the present case in view of the 
small pp scattering lengths, whose real parts are typically 
in the order of only one fermi. Our suspicion is nourished, 



in particular, by the observation that the experimental 
pp mass spectrum of the comparable decay J/ 1 fy—>ir pp 
does not show any obvious sign of a pp FSI, while ap- 
plicaton of the Watson-Migdal prescription would yield 
a similar, strong, near-threshold enhancement as for the 
7PP channel for any of the NN models we utilized. We 
demonstrated that a consistent qualitative description 
of the pp invariant mass spectra from the decay reac- 
tions J/^^^fpp and J/\I/— >7r°pp can be achieved, how- 
ever, within a more refined treatment of FSI effects as it 
follows from a DWBA approach. 

Though our study shows that the enhancement seen 
in the decay J/^^^pp could indeed be a result of the 
FSI in the pp system one has to admit that due to the 
uncertainties mentioned above and the controversal sit- 
uation in the ir°pp channel explanations other than final 
state effects cannot be ruled out at the present stage. 
Since the available data on the reaction J/^^ir^pp are 
afflicted by large error bars it would be desirable to ob- 
tain improved experimental information here that allows 
one to quantify the extent of pp FSI effects in this re- 
action. As already mentioned in the introduction, other 
reactions involving the pp system in the final state, such 
as B + ^K + pp decay [2] or B°^D°pp [3] do show some 
indications for pp FSI effects. But here too the quality 
of the available data is too poor to permit any reliable 
conclusions. 

In this context let us emphasize that pp mass spec- 
tra for the reactions J/^^upp or J/^^r/pp would be 
rather interesting for clarifying the role of pp FSI effects 
in the decay of the J/^> meson. Both reactions restrict 
the isospin in the pp system to be zero so that one could 
explore the FSI for specific pp partial waves, namely 11 So 
and 13 Si. 
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